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ABSTRACT 

High ambient interstellar pressure is suggested as a possible factor to explain the ubiquitous ob- 
served growth-rate discrepancy for supernova-driven superbubbles and stellar wind bubbles. Pressures 
of P/k ~ 10 5 cm -3 K are plausible for regions with high star formation rates, and these values are in- 
termediate between the estimated Galactic mid-plane pressure and those observed in starburst galaxies. 
High-pressure components also are commonly seen in Galactic ISM localizations. We demonstrate the 
sensitivity of shell growth to the ambient pressure, and suggest that superbubbles ultimately might serve 
as ISM barometers. 

Subject headings: galaxies: ISM — Magellanic Clouds — ISM: bubbles — ISM: general — supernova 
remnants 



1. INTRODUCTION 

Mechanical feedback from supernovae (SNe) and stel- 
lar winds is an important driver of galactic evolutionary 
processes. It affects the phase balance and physical con- 
ditions of the interstellar medium (ISM), which in turn 
determine star- forming conditions, galactic chemical evo- 
lution, and properties of the intergalactic medium. The 
standard paradigm for mechanical feedback is based on 
the model for adiabatic evolution of the shells and super- 
bubbles (e.g., Pikel'ner 1968; Weaver etal. 1977; Mac Low 
& McCray 1988) that are the direct consequence of SN 
and stellar wind action. 

While this model for pressure-driven superbubbles is 
broadly consistent with observations spanning scales from 
individual stellar wind bubbles to galactic superwinds, 
nagging problems persist in comparisons with observations 
(e.g., see Oey 2004 for a review). Specifics regarding the 
energy budgets, fate of the shock-heated 10 6 K gas, and 
later-stage evolution are lacking and have profound conse- 
quences for galactic evolution. 

One problem that is empirically well-established is 
the result that most superbubbles apparently grow more 
slowly than expected. This has been observed in individ- 
ual stellar wind bubbles such as Wolf-Rayet nebulae (Tref- 
fers & Chu 1982; Garci'a-Segura & Mac Low 1995; Drissen 
etal. 1995), as well as in superbubbles powered by OB as- 
sociations (e.g., Oey 1996a; Oey & Smedley 1998; Brown 
etal. 1995; Saken etal. 1992). This growth-rate discrep- 
ancy has been identified in young, nebular shell systems, 
in which the parent OB association is still present; thus 
the input mechanical power is well-constrained. The dis- 
crepancy is seen both in objects that show no evidence of 
previous supernova activity, and in ones where one or two 
supernovae have already exploded (Oey 1996a; hereafter 
096). 

For constant input mechanical power L and ambient 
number density n, the evolution of the shell radius is given 



by (e.g., Castor, McCray, & Weaver 1975), 



R = 68.9 (iss/n) 1 / 5 t 3 e /5 



pc 



(1) 



where L 38 is L in units of 10 38 ergs -1 , and t 6 is age of 
the bubble in Myr. The shell expansion velocity v is the 
time derivative of equation 1 . One possible solution to the 
growth-rate discrepancy suggests that the input parameter 
L/nis systematically overestimated. For eight nebular su- 
perbubbles with well-constrained R, v, L, and t, Oey and 
collaborators (096; Oey & Massey 1995; Oey & Smedley 
1998) showed that L/n would need to be reduced by a 
factor of several, perhaps up to an order of magnitude, to 
reconcile the observations with prediction. Since stellar 
wind power L is sensitive to the stellar mass, a substan- 
tial uncertainty in L is not unreasonable. As shown by 
multi-wavelength observations of three of the superbub- 
bles (Oey etal. 2002), the multi-phase ambient ISM also 
renders n similarly uncertain. However, the implication of 
a systematic growth-rate discrepancy remains difficult to 
explain. 

Another favorite candidate to solve the problem is cool- 
ing of the hot interior whose pressure drives the shell 
growth. If this scenario is correct, it implies a departure 
from the adiabatic evolution. The mass within the hot 
bubble interior is dominated by material evaporated from 
the cool shell walls, and could be supplemented by addi- 
tional material evaporated and ablated from small clouds 
that are overrun by the expanding shocks (Cowie & Mc- 
Kee 1977; McKee etal. 1984; Arthur & Henney 1996). 
This enhanced interior density would facilitate radiative 
cooling. Silich etal. (2001) and Silich & Oey (2002) also 
show that the enhanced metallicity caused by SN explo- 
sions and stellar products can further facilitate the cooling, 
especially for low- metallicity objects. However, increased 
X-ray luminosities that are expected from enhanced cool- 
ing thus far have not been observed (Chu et al. 2003; Chu 
etal. 1995). 

The superbubbles studied by Oey and collaborators are 
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all located in the Large Magellanic Cloud (LMC). For 
these objects, Silich & Franco (1999) suggested that the 
ambient environment and viewing geometry conspire to 
yield misleading observed shell dynamics. They suggest 
that the supcrbubblcs arc more extended perpendicular 
to the galaxy's plane, as would be expected in the plane- 
stratified density distribution of disk galaxies (but see also 
Maciejewski & Cox 1999). The elongation of the shells 
would not be apparent because of the LMC's almost face- 
on orientation. While this is an attractive suggestion for 
the LMC objects, it does not explain the growth-rate dis- 
crepancy seen in Galactic (e.g., Brown etal. 1995; Saken 
etal. 1992) and M33 (Hunter etal. 1995) objects. 

Nevertheless, it is apparent that the ambient environ- 
ment plays a crucial role in the superbubble growth and 
evolution. In addition to the work of Silich & Franco 
(1999), other studies have shown that the shell dynamics 
are sensitive to the ambient density structure (e.g., Oey 
& Smedley 1998; Mac Low etal. 1998). Multi-wavelength 
observations also show that the ambient multiphase gas 
distribution is difficult to constrain without direct such 
observations (Oey etal. 2002). 

2. THE INTERSTELLAR PRESSURE 

Continuing to focus on the ambient environment, this 
present work now investigates the effect of the ambient 
pressure. Since most superbubble growth is eventually 
expected to become confined by the ambient interstellar 
pressure (Oey & Clarke 1997), this parameter is also worth 
examining more closely. We note that Dopita ct al. (1981) 
suggested an active pressure mechanism to confine the 
shell of one LMC object by invoking exterior ram pressure 
caused by contraction of a surrounding interstellar cloud. 
Here, we suggest that the typical interstellar pressures in 
some systems may be higher than assumed. 

In our earlier work, we used a simple, semi-analytic, 
1-D model that integrates the shell's equations of motion 
(096; Oey & Massey 1995). For our sample of eight young, 
nebular superbubbles, we tailored the model input param- 
eters (L, n, t) according to the individual, empirically- 
derived values. As mentioned above, these highly con- 
strained models confirm the ubiquitous growth-rate dis- 
crepancy between the observed and predicted shell radius 
and expansion velocity (e.g., 096). The model predicts 
that the interior pressure evolves as (e.g., Weaver et al. 
1977), 

P l /k = 1.83 x 10 5 L^V/ 5 ^ 475 cm" 3 K . (2) 

In our above studies, the ambient interstellar pressure P e 
was estimated as P e = pc 2 s /"f. The soundspeed c s was usu- 
ally estimated as 10 kms -1 for ionized nebular gas, pro- 
vided that R remains smaller than the Stromgren radius; 
p is the mass density; and 7 is the equation of state index, 
which was taken to be 5/3 for the adiabatic condition. For 
n ranging between 1 and 10 cm~ 3 , as estimated for our 
objects, this yields P e /k = 9 x 10 3 to 9 x 10 4 cur 3 K. 

For the Milky Way, the total mid-plane ISM pressure is 
generally estimated to be around P/k ~ 3 x 10 4 cm~ 3 K. 
This includes roughly equal empirical contributions from 
the diffuse thermal pressure, magnetic field pressure, non- 
thermal velocity field or turbulent pressure, and cosmic 



ray pressure; a good discussion is presented by Slavin 
& Cox (1993). This value is consistent with the con- 
straint derived by Boulares & Cox (1990) that the pres- 
sure required to support the weight of the Galactic ISM is 
P/k ~ 2.8 x 10 4 cm" 3 K. 

In recent years, turbulent velocity structure is becoming 
more quantified as one of the major, and perhaps dom- 
inant, kinematic properties of the ISM in star- forming 
galaxies. This is linked to a new paradigm shift for the 
ISM to less distinct thermal phases, in which cool clouds 
are transient, unconfined features, rather than distinct 
and well-defined entities (Vazquez-Semadeni 2002; Krit- 
suk & Norman 2002). This turbulence-dominated view 
of the ISM implies a strong presence of components and 
localizations that are not in pressure equilibrium (Mac 
Low etal. 2001; Kim etal. 2001). Such components are 
now being ubiquitously identified in thermal pressure dis- 
tributions determined for lines of sight to Galactic stars 
(Jenkins & Tripp 2001; Wallerstein etal. 1995), where 
these are found to have thermal Pth/k ;> 10 5 cm -3 K. 
Also note that the turbulent pressure P t b/k ~ po^jk ~ 
7.5 x 10 4 cm -3 K for the hot, ionized medium (HIM) if 
we consider values of n ~ 0.05 cm -3 and turbulent veloc- 
ity dispersion a v ~ c s ~ 100 kms -1 . The magnetic field 
pressure also shows occasional hints for high-value com- 
ponents (e.g., Edgar & Cox 1993). On the other hand, 
the cosmic ray pressure contribution may be less relevant 
to confining superbubble growth since the shells may be 
transparent to cosmic rays (Slavin & Cox 1993). 

The above factors apply to the Milky Way, which is a 
giant disk galaxy whose midplane ISM pressure should be 
substantially higher than in a Magellanic irregular galaxy 
like the LMC, where all of our sample objects are located. 
Thus, one could argue that P e estimated by our earlier 
models are already on the high side of what might be ex- 
pected. On the other hand, several of the pressure terms, 
namely, the thermal, turbulent, and cosmic ray pressures, 
should scale with star formation rate (SFR) per unit vol- 
ume. Thermal pressures within star forming regions them- 
selves are of order P th /k ~ 10 5 - 10 6 cm -3 K (e.g., Mal- 
hotra etal. 2001), and superbubbles generally originate 
within such regions, although presumably they outgrow 
and outlive them. For a high filling factor or interstellar 
porosity generated by the superbubble activity, it is also 
quite likely that the magnetic pressure is also determined 
by the SFR. The magnetic field strength is correlated with 
gas density, yielding localized values of order a few hun- 
dred ^iG in Galactic molecular clouds and star-forming 
regions (Crutcher 1999). The expansion of superbubbles 
themselves also compresses the magnetic field, causing a 
self-induced impedance (Slavin & Cox 1992). On large 
scales, radio synchrotron measurements for other galax- 
ies are showing magnetic field strengths of ~ 15 ^G in 
actively star-forming galaxies, and up to 40 pG in spiral 
arms (Beck 2004). Such values imply magnetic pressures 
alone of order Ps/k ~ 10 5 cm -3 K. 

The LMC has one of the highest SFR per unit volume 
in the Local Group. Its interstellar porosity Q is exceeded 
only by that of IC 10 (Oey et al. 2001), although Q for the 
Milky Way is difficult to determine. Taking the SFR con- 
sistently estimated from the H II region luminosity func- 
tion for both the Galaxy (Oey & Clarke 1997) and the 
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Table 1 

lmc superbubble parameters 51 
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43 
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2 


2 
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49.342 
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3 


1 


7 


L301 


53 
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50.310 





(1) 


3 


3 


1 



a Data compiled by 096). Columns 6-11 represent numbers of stars in each mass bin; expected lifetime in Myr is 
shown in the column heading. 

Uncertainty ~ 10 - 15%. 

c Objects with ":" uncertain to 50%, but see text; others ~ 20%. Sec 096 for source references for v. 
d Uncertainty of order a factor of 2. 

°Valucs in parentheses show original number of stars implied by the IMF, from 096. 



LMC (Oey etal. 2001), we obtain porosities of Q ~ 0.3 
and 1, respectively. The formal Q for IC 10 is an order 
of magnitude higher yet. Thus it is not unreasonable that 
the LMC may have total interstellar pressures that are 
similar to, or in excess of, those of the Milky Way. We 
can also compare with the central few kpc of nuclear star- 
burst galaxies like M82, which have empirically determined 
P/k ~ 10 6 cm" 3 K (e.g., Lord etal. 1996). 'if the Milky 
Way ISM, with P/k ~ 10 4 cm -3 K, is typical for ordinary 
star-forming galaxies, then galaxies with SFR intermediate 
between "ordinary" and starburst, like the LMC, should 
be expected to likewise have interstellar pressures that are 
intermediate, namely, of order P/k ~ 10 5 cm- 3 K. Equa- 
tion 2 shows that such values would be significant in coun- 
teracting the pressure-driven shell growth. 

3. MODELS OF LMC SUPERBUBBLE EVOLUTION 

In light of the above, it is worth further exploring the 
effects of the ambient pressure parameter on the super- 
bubble evolution. We return to the same LMC objects 
studied by 096, now modeling these with 1-D hydrody- 
namical computations using the magneto-hydrodynamic 
fluid solver ZEUS-3D (version 3.4). This code is the up- 
dated, 3-D version of the two-dimensional code ZEUS-2D 
(Stone & Norman 1992), an Eulerian explicit code which 
integrates the equations of hydrodynamics for a magne- 
tized ideal gas. The code also works efficiently in one 
dimension, and we perform the simulations in spherical 
coordinates, with 1000 zones in the radial direction r. Tak- 
ing the time-dependent stellar wind mechanical luminosity 
L(t) = \ M(t) for each object from 096, we assume 
a constant terminal velocity Vqq = 3000 kms -1 , yielding 
a wind density p(t) — M(i)/47rr 2 i>oo, where M(t) is the 
wind mass-loss rate. 096 estimated L{t) from the observed 



and inferred stellar population (Table 1, below) based on 
empirical relations between wind properties and spectral 
type, and evolutionary models. The time-dependent wind 
conditions are set within the first few radial computational 
zones, centered at the origin. 

ZEUS-3D does not include radiation transfer, but we 
have implemented a simple approximation to derive the 
location of the ionization front for arbitrary density dis- 
tributions (see Bodenheimer et al. 1979). This is done by 
assuming that ionization equilibrium holds at all times, 
and that the gas is fully ionized inside the HII region. The 
position of the ionization front is given by J n 2 (r)r 2 dr « 
Q° /Anas, where Q° is the stellar H- ionizing emission rate 
and as is the Case B recombination coefficient. 

The models include the Raymond & Smith (1977) cool- 
ing curve above 10 4 K. For temperatures below 10 4 K, the 
shocked gas region is allowed to cool down with the radia- 
tive cooling curves given by Dalgarno & McCray (1972) 
and MacDonald & Bailey (1981). Finally, the photoion- 
ized gas is always kept at 10 4 K , so no cooling curve is 
applied to the HII regions (unless there is a shock inside 
the photoionized region). 

We assume an initially homogeneous ISM, which sup- 
plies the external pressure. Given the existence of large 
non-thermal "turbulent" velocities, of several km s _1 , cos- 
mic rays and magnetic fields, the total pressure should 
be obviously a combination of all of them, as discussed 
above. However, to simplify computation, we model the 
total ambient pressure as originating entirely from ther- 
mal pressure. Thus, in order to have thermal pressures 
of P/k = 1 x 10 5 cm -3 K, we adopt for the ISM density 
n <~ 16.7 cm- 3 and 6000 K for the temperature. This gives 
us a sound speed of c s .isoth ~ 7 kms . For a second set 
of models with P/k = 1 x 10 4 cm -3 K, we use an ISM 
density of n ~ 1.67 cm -3 and 6000 K for the temperature. 
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Fig. 1. — Modeled density profiles for the six LMC superbubbles, assuming an ambient P/k = 1 X 10 5 cm 3 K. The input stellar popula- 
tions are given in Table 1, along with observed parameters. The observable nebular shell is delineated by the vertical dashed and dotted line. 




Fig. 2. — Predicted velocity profiles for the models shown in Figure 1, with lines as in Figure 1. 



6 



Ambient Pressure and Superbubble Evolution 



Table 1 gives the observed nebular radius R, nebular 
expansion velocity v, and Q° in columns 2-4, respec- 
tively, for the sample objects. Columns 5-10 give the ob- 
served and pre-SN inferred massive star population down 
to 20 Mq. 096 describes how the input mechanical power 
L(t) due to stellar winds and SNe is estimated from the 
individual massive star populations as a function of time. 
Figures 1 and 2 of 096 show that L(t) for these objects 
typically have early values of 10 37 — 10 39 ergs -1 . All 
parameters are listed in Table 1 as measured and com- 
piled from the literature by 096, except for Q°, which is 
estimated by Oey & Kennicutt (1997). The objects are 
divide into "Pre-SN" and "Post-SN" categories: the Pre- 
SN objects show no enhanced X-ray emission or enhanced 
[S II] /Ha ratios, thus are presumed to have no prior SN 
activity. Objects listed as Post-SN show both enhanced 
X-ray emission (Chu & Mac Low 1990; Wang & Helfand 
1991) and [S IlJ/Ha, indicating recent SN impacts on the 
shell walls; the masses of the SN progenitors are estimated 
from the stellar initial mass functions (IMF) and included 
in L(t) for the models, as done by 096. 

Figure 1 shows our new models for the gas mass density 
as a function of radius, using the input L(t) determined by 
096 from the observed stellar population (Table 1) of the 
individual LMC objects. These models adopt a high ambi- 
ent pressure P/k = 1 x 10 5 cm~ 3 K. As discussed earlier, 
contributions from turbulent, magnetic, and multi-phase 
thermal pressure terms would presumably contribute to 
such high ambient values. The models are those corre- 
sponding to the observed radii in Table 1, with Pre-SN 
objects shown in the left column, and Post-SN objects 
shown on the right. Figure 2 shows the gas velocity as a 
function of radius for the same models. The location of 
the observable ionized nebula is shown between the verti- 
cal dashed and dotted line in both Figures. 

Comparison with Table 1 shows that the high-pressure 
models are in good agreement with the observed param- 
eters in all cases: the predicted ionized radius, expansion 
velocity, and age are all reasonably consistent with the 
observed values. In Figure 4 we reproduce Ha images of 
the supcrbubbles from Figure 1 of Oey (19966). Figure 4 
also shows that the putative Pre-SN objects DEM L106 
and L226 show thicker nebular shells, of order 10% of R, 
whereas the post-SN objects show compressed, filamen- 
tary morphology. DEM L31, presumed to be Pre-SN, is 
also quite filamentary, and the model also reproduces this 
structure, owing primarily to its high input power (Ta- 
ble 1). This is fully consistent with the predicted mor- 
phology of the ionized regions in the models. 

The velocity profiles shown in Figure 2 also imply ob- 
served velocities that are reasonably consistent with the 
data (Table 1). The latter are especially complex for the 
most filamentary objects, showing velocity components 
that vary by factors of 2-3. This is consistent with the 
extreme gradients in both density and velocity seen in the 
models. Figure 3 shows the resulting predicted, positive 
velocity profiles with respect to the systemic velocities, in- 
tegrated through the centers of the superbubbles. DEM 
L31 and DEM L301 show a significant range in velocity, 
and maximum values consistent with those observed. Pre- 
dicted velocities for the remaining objects are also consis- 
tent with the observed upper limits. We note that both ob- 
servations and models are likely to be sensitive to details of 



geometry and ionization in the post-SN objects, since com- 
plex shock structures are generated on short timescales. 
Thus it is unsurprising that the model velocity profiles 
tend to be more simplistic than the observations. We also 
see that the models for DEM L25 and DEM L50 imply that 
the ionization front does not penetrate the high-density 
shell (Figure 1); however, observations of these objects 
do show high-density, filamentary structure that is fully 
ionized for DEM L25, and mostly ionized for DEM L50 
(Oey et al. 2002). Oey & Kennicutt (1997) find excess 
in the nebular emission, by factor of a few, beyond what 
can be attributed to stellar photoionization. Thus, shock 
excitation is likely to enhance the ionized mass and radius 
of these objects. Nevertheless, the broad agreement in the 
dominant shell parameters across the sample demonstrates 
that a high ambient pressure alone can solve the growth- 
rate discrepancy that is widely observed in such objects. 

Figure 5 shows the density profile evolution modeled for 
the parameters of DEM L226 and DEM L25. The former 
is representative of Pre-SN, low-L objects; and the latter 
is representative of Post-SN, higher- L objects. The curves 
show a sequence of density profiles over 10 Myr at 4 x 10 5 yr 
intervals, from top to bottom. We see that the evolution of 
both objects is similar, but that the model for DEM L226 
does not generate enough internal pressure to maintain 
shell growth within the modeled time frame: the central 
cavity is seen to become pressure-confined, and to start 
collapsing. On the other hand, the predicted evolution for 
DEM L25 does maintain shell growth within the modeled 
time frame. In both models, the shell thickens and diffuses 
as the expansion falls below the soundspeed, thereby allow- 
ing the shock front and shell to dissipate. This is consistent 
with behavior found by Slavin & Cox (1992) and Garcfa- 
Segura & Franco (1996) for similar pressure-confined mod- 
els. At similar, late, subsonic/pressure-confined stages, 
the structure is vulnerable to disruption by ambient turbu- 
lence, especially if the latter is an important contributor to 
a high-pressure situation. In the meantime, both models in 
Figure 5 also clearly show compression in the shell caused 
by supersonic SN impacts around 4.1 Myr. Additional, 
subsequent SN impacts also can be seen in the models. 
We do note that DEM L25 is the only object whose age 
in Figure 1 is somewhat inconsistent with the predicted 
stellar population. For a standard IMF, one 40 M star is 
expected to have existed, which presumably has since ex- 
ploded, implying an age of 5.3 - 7.8 Myr. The evolutionary 
sequence in Figure 5 shows R ~ 53 pc, vs the observed 43 
pc, at that stage. This is not a large discrepancy, and the 
existence of the putative 40 M Q SN progenitor is also ex- 
tremely uncertain in view of the stochastic effects in the 
IMF. 

We also compare with similar models assuming a lower 
ambient pressure of P/k = 1 x 10 4 cm~ 3 K, implemented 
as described above. Figures 6 and 7 show the models cor- 
responding to the observed radii for the same six objects. 
While this ambient P/k value is more consistent with stan- 
dard expectations, it is apparent that the models are in 
poor agreement with the data. The collective ages are sys- 
tematically too young, thus demonstrating the ubiquitous 
growth-rate discrepancy, and the predicted shell parame- 
ters are more difficult to reconcile with the observations. 
Since none of the objects have reached the SN stage, the 
Post-SN objects retain thick, extended shells. Indeed, all 



M. S. Oey & G. Garci'a-Segura 



600 



400 



~| 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

D31 

2.6 Myr 



S 200 





800 
600 



o 

o 400 

Oh 




2 



200 



400 



o 
ft 



200 



I 1 1 1 1 1 1 1 1 1 h 



D106 

2 Myr 



I 1 1 1 1 1 1 1 1 1 h 



D226 
2.8 Myr 




i i i i i i i I i i i i i i i i i i i i i i i i i i i i i i 



~i 1 1 1 1 1 1 1 1 r 



-| 1 ' 1-=1 25 



H 1 1 1 1 h 



\ 1 1 1 1 1 1 1 1 1 h 



D301 
4.6 Myr 




D25 

4.4 Myr 



20 
15 
10 



H 1 1 1 1 h 

D50 

5.6 Myr 



10 



1000 



500 



_L 



_L 



20 40 

velocity (km/s) 



60 



20 40 60 

velocity (km/s) 



Fig. 3. — Predicted positive Ha line profiles for the models shown in Figure 1. Note that the total line profile will be symmetric around 
the origin (systemic velocity). 
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Fig. 4. — Ho images of the LMC superbubbles, from Oey (1996b). Each image is 16/67 square. 




Fig. 5. — Predicted evolution of the density profiles for DEM L226 and DEM L25. Modeled radial profiles arc shown at intervals of 4 X 10 5 
yr, over a total of 10 Myr. The observed R are consistent with the models at ages of 2.8 Myr and 4.4 Myr for DEM L226 and DEM L25, 
respectively. 
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Fig. 7. — Predicted velocity profiles for the models shown in Figure 6. 
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Fig. 8. — Predicted positive Ha line profiles for the models shown in Figure 6. 
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but one of the objects are predicted to have a surrounding 
ionized halo that is at most a factor of 3 lower in density 
than the inner shell; these halos arc not observed. An- 
other especially interesting feature in these models is the 
appearance of a double shell structure: a second, neutral 
shell is seen at roughly 2-3 times the radius of the ionized 
inner shell. This outer shell results from the expansion of 
the photoionized H II region in a D-type ionization front. 
Oey et al. (2002) presented H I observations of the three 
post-SN objects (DEM L25, DEM L50, and DEM L301), 
at a spatial resolution of 50" (12.5 pc); they found no ap- 
parent evidence of such secondary shells. The existence of 
neutral secondary shells could potentially offer a diagnos- 
tic of lower ambient pressures. Figures 7 and 8 also show 
the predicted velocity structure and Ha line profiles for 
these young objects, which, as expected, are less consis- 
tent with the observations in Table 1, than are Figures 2 
and 3. 

4. CONCLUSION 

Our models clearly show that increasing the ambient in- 
terstellar pressure by an order of magnitude, from P e /k = 
1 x 10 4 to 1 x 10 5 cm~ 3 K, can impede the shell growth to 
a degree that could fully account for the observed growth- 
rate discrepancy. In §2, we presented arguments that such 
high interstellar pressures could exist, especially based on 
the dependence of P e on star-formation rate. While other 



factors mentioned in §1, namely, overestimated L/n, ele- 
vated radiative cooling, and viewing geometry, could all be 
additional factors that contribute to the growth-rate dis- 
crepancy, we note that the multi-phase gas morphology is 
more consistent with high interstellar pressure dominating 
this effect. 

Finally, as noted by Oey & Clarke (1997), the assumed 
global value of P e plays a critical role in determining the 
characteristic final sizes of old, SN-dominated superbub- 
bles, and hence, the superbubble size distribution, which is 
dominated by pressure-confined shells. This, in turn, de- 
termines the interstellar porosity and filling factor of the 
hot, ionized medium in star-forming galaxies. With ad- 
equate clarification in the superbubble evolution process 
and input parameters, the superbubble sizes, kinematics, 
and morphologies could potentially provide barometers for 
the interstellar pressure. These diagnostics could be espe- 
cially useful in other galaxies, which have fewer available 
empirical pressure indicators than the Milky Way. 

We are grateful to the anonymous referee, whose com- 
ments led to a much stronger analysis of this problem. 
Thanks also to Dave Strickland for comments on the 
manuscript and for prodding us to finally carry out this 
work. We are also pleased to acknowledge comments and 
discussions with Don Cox and Robin Shelton. Much of this 
work was carried out by the PI at Lowell Observatory. 
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